Generic gauge fields in the Hubbard model: emergence of pairing interaction 
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The spin-rotationally invariant SU(2) approach to the Hubbard model is extended to accom- 
modate the charge degrees of freedom. Both U(l) and SU(2) gauge transformation are useed to 
factorize the charge and spin contribution to the original electron operator in terms of the emergent 
gauge fields. It is shown that these fields play a similar role as phonons in the BCS theory: they 
provide the "glue" for fermion pairing. By tracing out gauge bosons the form of paired states is 
established and the role of antiferromagnetic correlations is explicated. 
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To understand the physics of strongly correlated (SC) 
systems one frequently employs a slave particle (SP) de- 
composition of the electron operator ji£ Analytic theory 
behind this formulation hinges on the treatment of the 
constraint of no double occupation based on the assump- 
tion that the on-site interaction energy U can be renor- 
malized to infinity. The redundancy in representations 
used to enforce the constraint naturally leads to vari- 
ous gauge theories £ It is sometimes supposed that slave 
particles can be liberated at low energies, in which case 
the slave-boson and fermion degrees of freedom take on 
a physical meaning with the spin-charge separation as a 
result. However, it was shown that gauge theories as- 
sociated with SP representations of correlated electrons, 
such as the t — J model, are always confining. 4 The rea- 
son is that the slave-particle gauge theory is infinitely 
strongly coupled- there is no intrinsic kinetic energy for 
the gauge field. Therefore, the slave particle represen- 
tations of correlated electron models are not belonging 
to the class of generic gauge theories. Although mat- 
ter fields which couple to the non-generic SP gauge field 
might generate gauge dynamics in the low-energy effec- 
tive action^ the problem remains: how to derive emer- 
gent gauge fields from the microscopic formulation of the 
theory of SC electrons and not just to write down a La- 
grangian that contains them - as in the SP approach. 
This issue is of basic concern because the understanding 
of the mechanism of superconductivity in cuprates re- 
quires the knowledge of bosons mediating the pairing as 
well as the nature of the paired states. Here, the underly- 
ing attraction force appears very puzzling since it is hard 
to reconcile the microscopic attractive interaction with 
the completely repulsive bare electron-electron forces. 

In the present paper we extend the SU(2) spin- 
rotationally invariant approach to the Hubbard model, 6 
which makes no assumptions regarding the magnitude of 
the Coulomb energy U, to accommodate on equal foot- 
ing the charge degrees of freedom. Using U(l) and SU(2) 
transformation we explicitly factorize the charge and spin 
contribution to the original electron operator in terms of 
the corresponding emergent gauge fields. We show that 
these fields play a similar role as phonons in the BCS the- 
ory: they provide the "glue" for fermion pairing in the 
SC system. By tracing out gauge bosons we explicitly 



calculate and the form of paired states and explicate the 
role of antiferromagnetic (AF) correlations. 

Our starting point is the purely fermionic Hubbard 
Hamiltonian Tt = Tit + T~iu'- 



H = -tJ2 [4(r)c Q (r')+h.c.]+5]C/n T (rK(r). (1) 
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Here, (r, r') runs over the nearest-neighbor (n.n.) sites, 
t is the hopping amplitude, U stands for the Coulomb 
repulsion, while the operator cj,(r) creates an electron 
with spin a =|,| at the lattice site r. Furthermore, 
n(r) = riffr) + n^(r) is the number operator, where 
n Q (r) = c^(r)c a (r). Usually, working in the grand 
canonical ensemble a term — fJ-^2 r n(r) is added to TL in 
Eq. Q with [l being the chemical potential . It is custom- 
ary to introduce Grassmann fields, c«(rr) depending on 
the "imaginary time" < r < (3 = l/k B T, (with T being 
the temperature) that satisfy the anti-periodic condition 
c a (rr) — — c Q (rr + /3), to write the path integral for the 
statistical sum Z = J [DcDc] e~ 5 I c ' c l with the fermionic 
action 

S[c,c]=S B [c,c] + / dTH[c,c\, (2) 
Jo 

that contains the fermionic Berry term 
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B[C,C\ = 
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drc a (rr)d T c a (rr) . 



(3) 



For the SC system it is crucial to construct a covariant 
formulation of the theory which naturally preserves the 
spin-rotational symmetry present in the Hubbard Hamil- 
tonian. For this purpose the density-density product in 
Eq.Q) we write, following Rcf 6, in a spin-rotational in- 
variant way: 



Hi 
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2 (rr) - [fi(rr) • S(rr)]' 



(4) 



where S°(rr) < = 5E aa '4HC' c «'( r y) denotes 
the vector spin operator (a = x, y, z) with a a be- 
ing the Pauli matrices. The unit vector f2(rr) = 



2 



[sin#(rr) cos y(rr), sin-#(rT) sin y(rr), cos $(rr)] written 
in terms of polar angles labels varying in space-time spin 
quantization axis. The spin-rotation invariance is made 
explicit by performing the angular integration over f2(rr) 
at each site and time. By decoupling spin and charge 
density terms in Eq.(@J) using auxiliary fields q(tt) and 
iV(rr) respectively, we write down the partition function 
in the form 



Z = J [Vn] J [VVVq] J [DcDc] x 



x e 



-s[ny,g,c,c] 
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where[2?ft] = ]J rTk ^(^ k )d^rr k )d v (rr k ) is the spin _ 
angular integration measure. The effective action reads: 



rP 

S[n,V,g,c,c] = V / dr 
„ Jo 



q 2 {tt) V 2 {vt) 
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+ iV{rT)n(YT) + 2q(yt)VI(yt) ■ S(rr)] 
f P 

+ S B [c,c}+ / dTH t [c,c\. (6) 
Jo 

Simple Hartree-Fock (HF) theory will not work for a 
Hubbard model in which U is the largest energy in 
the problem. One has to isolate strongly fluctuating 
modes generated by the Hubbard term according to the 
charge U(l) and spin SU(2) symmetries. To this end 
we write the fluctuating "imaginary chemical potential" 
iV(rr) as a sum of a static Vo(r) and periodic function 
V(rr) = Vb(r) + V(rr) using Fourier series 



V(rr) 



OO 
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with u„ = 27rn//3 (n = 0, ±1,±2) being the (Bose) Mat- 
subara frequencies. Now, we introduce the U(l) phase 
field 4>{tt) via the Faraday-type relation 



<90(rr) 
dr 



V(tt). 



(8) 



Since the homotopy group tti[J7(1)] forms a set of inte- 
gers, discrete configurations of </>(rr) matter, for which 
0(r/3) - </>(r0) = 27rm(r), where m(r) = 0,±1,±2, ... 
Thus the decomposition of the charge field V(rr) con- 
forms with the basic m = topological sector since 
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field, except for the zero frequency. Subsequent SU(2) 
transformation from f a (rr) to h a (rr) operators, 



(10) 



with the constraint |(i(rr)| 2 + |C2( rr )| 2 = 1 takes away 
the rotational dependence on ri(rr) in the spin sector. 
This is done by means of the Hopf map R(rr)o' z R'l'(rT) = 
<x ■ fi(rr) that is based on the enlargement from two- 
sphere 5*2 to the three-sphere S3 ~ SU(2). The unimod- 
ular constraint can be resolved by using the parametriza- 
tion 



Ci(rr) = e ~5[¥ , ( rr )+x( I - r )] 



cos 



tf(rr) 



C 2 (rr) 



= e iM r T)-x(rr)] gin 
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with the Euler angular variables ip(rr), $(rr) and x(rr), 
respectively. Here, the extra variable x( rr ) represents 
the U(l) gauge freedom of the theory as a consequence 
of S2 — > 53 mapping. One can summarize Eqs JHJ and 
HI U| l by the single joint gauge transformation exhibiting 
electron operator factorization 



c Q (rr) = ^z{YT)R aa ,(rT)h a >(YT) 



(12) 



where R(rr) 



e -idzip{rT)/2 e - 



■ia y d(rT)/2 p -ia zX (rT)/2 



unitary matrix which rotates the spin-quantization axis 
at site r and time r. Eq. l|12H reflects the composite nature 
of the interacting electron formed from bosonic spinorial 
and charge degrees of freedom given by R aa '(rr) and 
z(rr), respectively as well as remaining fermionic part 
h a (rr) . In the new variables the action in Eq. © assumes 
the form 



fP 

S [fl, (j>, g, h, h] = Sb[}i, h] + / dTHn.<f>[p,h,h] 
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+S [0] + 2 V / drA(rr) • S h (rr), (13) 
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tffarr) = JqV{vt) = 0. Furthermore, by perform- where S ^( rr ) = 2 E Q7 Mrr)<x Q7 /i 7 (rT). Furthermore, 



ing the local gauge transformation to the new fermionic 
variables / Q (rr): 



c Q (rr) 




_ c Q (rr) _ 





z(rr) 
z(tt) 



fa(TT) 

/a(rr) 
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where the unimodular parameter |2(rr)| 2 = 1 satis- 
fies z(rr) = e l ^ rT \ we remove the imaginary term 
i Jq drV (rr)n(rT) for all the Fourier modes of the V(rr) 



stands for the kinetic and Berry term of the U(l) phase 
field in the charge sector. The SU(2) gauge transforma- 
tion in Ea. Hl(J|) and the fermionic Berry term in Eq.Q 
generate SU(2) potentials given by R^(rr)9 r R(rr) = 



3 



rr ■ A(rr), where 
A x {vt) = 



1 . 

-i?(rr) sin x(rr) 

1 

—<p(rr) sin 9(rr) cos x( rr ) 



2 . 

A y (vr) = — 7?(rr) cos x( rr ) 

+ -(yj(rr) sin 6>(rr) sin x(rr) 

z z 



(15) 



In analogy to the charge U(l) field the SU(2) spin system 
exhibits emergent dynamics. By integrating out fermions 
the last term in Eq. I|13|) will generate the kinetic term for 

the SU(2) rotors 5 [«] = -(l/£ s ) E r la ^A(rr)- A(rr) 
in the form 



So[n] 



i 



eh 



i? 2 (rr) +(^ 2 (rr) 



+ x 2 ( rr ) + 2<^(rr)x(rr) cos i?(rr)l 
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with £ s being of order of U close to half filling. The 
first order term in A(rr) fields gives rise to the usual 
spin Berry contribution. If we work in Dirac "north 
pole" gauge x( rr ) = — L f( rT ) one recoverss the famil- 
iar form S B [«] = -i(p/U)J2 r So ' *-<p(rr)[l - costf(rr)]. 
The fermionic sector, in turn, is governed by the effective 
Hamiltonian 

rin^ = ^£»(rr)[/i T (rT)/i T (rr) - ^(rr)/i;(rr)] 

r 

-t ^ z{vt)z{v't) [Rt(rr)R(r'r)] ^(rr)/i 7 (r'r) 



(r,r'> 



rex 



(17) 



where /2 = /x — nfU/2 is the chemical potential with 
a Hartree shift originating from the saddle-point value 
of the static variable Vb(r) with = nu\ + ^fej, and 
ri/m = (ji a (vT)h a {vT)) . The chief merit of the gauge 
transformation in Ea. (|12f> is that we have managed to 
cast the SC problem into a system of non-interacting h 
fermions submerged in the bath of strongly fluctuating 
U(l) and SU(2) gauge potentials coupled to fermions via 
hopping term plus Zeeman-type contribution with the 
massive field q(tt). In the AF phase, at the half- filling, 
it assumes the staggered form q(tt) = A c (— l) r with A c 
being the charge gap A c ~ U/2 for U/t ^> 1. However, 
a nonzero value of A c does not imply the existence of 
AF long-range order. For this the angular degrees of 
freedom O(rr) have also to be ordered, whose low-lying 
excitations are in the form of spin waves. 

It is well known that phonons play the role of the 
"glue" that is responsible for the formation of Cooper 
pairs in conventional superconductors. Now we show 
that U(l) and SU(2) emergent gauge fields, the collec- 
tive high energy modes in the SC system, take over the 
task which was carried out by phonons in BCS super- 
conductors. In a way similar to phonons these gauge 
fields couple to the fermion density type term via the 
amplitude t, see Ea. l|17[) . Now we evaluate the effective 
interaction between fermions by tracing out the gauge 
degrees of freedom. To this end we write the partition 
function as Z = f[DhVh]e- s[h ' h \ where 



S[h,h] = -In J [VSVD(f>], 



-S[ii,<p,h,h] 



(18) 



generates cumulant expansion for the effective fermionic 
action. We concentrate on the second order term in the 
hoping amplitude t containing four fermion operators: 



\ |ri — |— n.n. oca.' 

x ]T z(r 2 r0<r^r0E[ Rt ( r ^0R(4r0] 77 ,ft 7 (r 2 r0V(4r0y 

|r 2 — r' 2 \=n.n. 77' / 



(19) 
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where (...) denotes averaging over U(l) and SU(2) gauge 
fields. The averaging in the charge sector is performed 



with the use of the U(l) phase action in Ea. (|14|l to give 
(2(r 1 T)z(r / 1 r)2(r 2 T y )z(^r')) 

- (S riyi S r2 y 2 + S ri y 2 S r ' i . r2 ) X 



u 



(20) 



Equation (|20|1 reflects the local (in space) nature of 
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charge excitation and contains only the non-topological 
part of the four-point charge correlator. Away from 
half-filling the dynamics of spin variables is slower as 
compared to the charge counterparts, allowing to treat 
SU(2) variables as local in time R(rr') = R(rr) + (r' — 
r)5 r R(rr)+0[(r'-T) 2 ]. Furthermore, in the low temper- 
ature limit (on the energy scale given by U), by making 
use of the formula 



lira 

T->0 



we arrive at 



dr'e 



v-t'\U 
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(21) 



SW[h,h] = ~ f dr V y)[fc a (r'r)V(rr)] t 

U J ° <r,r<> <W 

x (A/ Q / Q;77 /(rr;r'r|r'TrT)) /i 7 (rV)/iy (rr), (22) 

where (...) denotes averaging over the remaining spin- 
angular variables and 



[Rt(rr)R(r'r)] a , a [Rt(r'r)R(rr)] (23) 



Now, employing the composition formula for rotational 
matrices^ 

R t (rr)R(r'r) 



1 

71 



e f*T_ 



(24) 



with T±(rr,r'r) = ^1 ± fi(rr) • n(r'r), where $ = 
$[r2(rr), fi(r'r), z] is the signed solid angle spanned 
by the vectors O(rr), O(r'r) and z with $ = 
<&[Sl(rr), — fi(r'r)] — 2y?(rr), we finally conclude that 



5' 2 '[fc,ft] = -V / dr [j^(rrr'T)^(rrr'T) 



where 



and 



J^JFCrrrV^rrrY)] , 



2t 2 

JA/ F = ir [lT<n(n-)-n(rV)>] 



(25) 



(26) 



■A(r-rrV) = 



h,(rT)h l (r'T)-h l (rT)hi(r'T) 
V2 



, . . ft t (rr)/i T (r'r) + fti (rr)/ii (r'r) , N 

^(rrr'r) = n ; n ; U ; u - (27) 

v2 

being the valence bond operators^ The rotational in- 
variance of the right-hand side in Ea. (|25|) is manifest 
since — A(yty't)A(yty't) = S/,(rr) ■ S/j(r'r) — \ and 
T (yty' t)T {yty' t) — S/j(rr) • S/^r'r) + |, respectively. 
The effective non-retarded interaction J a > in front 
of the ^4(rrr'r)^l(rrr / r) term constitutes the attractive 
potential for fermion pairing. By HF decoupling of the 



four-fermion term in Eq. (|25|l one obtains the "ci- wave" so- 
lution for the singlet pairing gap Ad(rr') = (A(yty't)}. 
The role of antiferromagnetic correlations is also appar- 
ent due to the presence of spin-angular correlation func- 
tion in Ea. H26|) . For example, in a fully developed AF 
background f2(rr) • n(r'r) = — 1 for r and r' on neigh- 
boring sites, so that J a = 4t 2 jXJ and Jp = 0, thus pro- 
moting "d-wave" pairing. The expectation value of the 
second valence bond operator {T(yty't)) competes with 
Arf(rr'), since it enhances the kinetic energy of fermions, 
which eventually results in suppression of the d-wave gap 
at higher doping level. The actual strength of the ef- 
fective interaction in Ea. H26|) requires that the quantity 
(f2(rr) • ri(r'r)) has to be determined self-consistently. 
To this end, in a way similar to Ea. (|18fl . one should in- 
tegrate from the effective Hamiltonian in Ea. i|17|) charge 
and fermion variables to obtain the action for the spin 
rotational degrees of freedom: 

S W = tE f dr[n{YT)-n{Y'r)-l] 

(rr') J ° 



(rr') 



S c (rr) • S c (r'r) - - 



(28) 



where we made use of the formula fl = itr(R<7 z Rt<x), 
while the variables (rr) = \ J2 a -y Ca(^T)o- ai ( 1 (YT) are 
the "bosonic" spins in the complex-projective (CP 1 ) for- 
mulation, see Ea. 1)11(1 . Here, J = -£t(«/iT — n hi) 2 in- 
dicates that Coulomb energy U induced Hubbard band 
splitting is a necessary prerequisite to sustain AF corre- 
lations. Although, the fermions h a (Yr), h a (Yr) play the 
role similar to "spinons" in the slave particle formula- 
tion a major quantitative difference appears. In SP the- 
ory at half-filling "spinons" are paired at T — with 
Ad 0. This is clearly impossible here due to the pres- 
ence Zeeman-type band splitting term in Eci. (|17|) that 
marks the onset of charge gap. It prevents a non-zero 
solution for since at half filling p is approaching the 
high energy charge gap, while A^ is governed by the ex- 
change energy J a <C U. We ascribe this discrepancy to 
the inherent inability in the SP scheme to give an ac- 
count of the high-energy effects that are the hallmark 
of the Mott physics & The reason is that the SP theory 
handels exclusively with with a strictly low- energy effec- 
tive Hamiltonian. This feature is odd e.g. with exper- 
iments which show that the superconducting transition 
in cuprates is accompanied by changes in the optical re- 
sponse, even at energies of the order of 100 times the 
critical temperature, clearly pointing out the importance 
of the high-energy effects on the scale given by Fi- 
nally we observe that superconductivity demands more 
than just paired fermions - it also requires phase coher- 
ence in the charge sector distinguished by the variables 
z(yt) — e l ^ TT \ Therefore a fully self-consistent theory 
requires a counterpart of the action in Ea. (|28|l for the 
phase variables and phase stiffnesses that are responsible 
for the actual superconducting state. 11 
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